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ABSTRACT 

A number  of  techniques  suggested  for  numerical  analysis  of  scattering  from  2D-objects 
with  complex  border  shape  and/or  arbitrarily  changing  electromagnetic  parameters  (for 
brevity,  from  2D-objects  of  arbitrary  composition)  is  very  large  indeed,  even  if  to  take 
into  account  only  the  approaches  based  on  finite  element  method  (MFE).  Nevertheless, 
some  methodological  and  “technological”  questions  of  MFE  application  have  been 
under  consideration  ever  since  first  implementations  of  the  method.  The  task  of  correct 
field  behavior  description  near  the  border,  where  the  electromagnetic  parameters  change 
abruptly,  belongs  to  the  first  group  of  questions.  The  convenient  means  of  problem 
statement,  such  as  specification  of  scattering  object  geometry  and  its  electromagnetic 
parameters  variation,  is  an  important  “technological”  task. 

The  aim  of  this  paper  is  to  present  the  MFE-based  numerical  procedure  to  evaluate  the 
RCS  signature  of  2D-objects  of  arbitrary  composition,  including  dielectric,  metal  and 
plasma-like  objects.  The  scattering  of  waves  of  two  polarizations  (H-wave  with  E. , 

Hr  components,  and  E-wave  with  H. , Er  components)  is  calculated.  The  original 

variants  of  solutions  to  the  two  above-mentioned  questions  were  implemented:  the 
weighted  residual  minimization  type  condition  is  implemented  on  the  border  of 
scattering  object  in  order  to  achieve  the  continuity  of  tangential  to  the  border  field 
components,  and  AutoCAD  is  used  for  data  input. 

CHARASCTERISTIC  FEATURES  OF  THE  NUMERICAL  PROCEDURE 


Wave  equations  for  E.  or  H components  are  solved  with  the  nodal  FE.  The  objects 

with  arbitrary  cross  sections  are  placed  inside  the  rectangular  region,  where  the  MFE  is 
used.  The  rectangularity  permits  a simple  adjustment  of  the  mesh  with  respect  to  the 
objects’  borders. 

The  AutoCAD  possesses  the  convenient  user’s  interface;  it  is  widely  used  in  industry. 
That  is  why  AutoCAD  was  chosen  for  inputting  data  on  the  geometry  of  the  problem 
and  electromagnetic  parameters  distribution.  The  borders  of  objects  are  outlined  with 


cubic  spline  r(u ) = a0  + dj  • u +a2 
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electromagnetic  parameters  can  either  be  set  analytically  or  as  a surface  z(x,  y)  above 
the  cross  section.  The  data  are  read  from  AutoCAD  output  file,  and  then  the  mesh  is 
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adjusted  with  the  use  of  parametrical  object's  border  representation,  see  the  examples  of 
AutoCAD  drawing  and  automatically  generated  mesh  in  Figl  a.b). 


| 


a)  b) 

Fig.l 

The  distribution  of  E.  or  H _ is  found  as  a solution  of  the  problem.  In  comparison  with 
the  known  edge-based  finite  element  technique  [1]  the  reduction  of  resulting  SLAT 
order  is  achieved,  while  the  process  of  SFAF  creation  is  simpler. 

The  question  of  the  tangential  field  component  (Hx  for  1 1-wave  or  /C  for  F-wavc)  on 

the  border  of  two  media  is  solved  by  the  corresponding  boundary  condition  imposition, 
analogously  to  [2],  but  again  with  respect  to  the  only  z-component  of  the  field. 
Generally  the  equation  of  SFAF  has  the  form  R\  +a  • Rj  = 0,  where  R\  and  /?2  arc 

the  weighted  residuals,  corresponding  to  satisfying  wave  equation  and  tangential  field 
component  continuity  on  the  border  of  two  media,  a is  a coefficient.  The  calculated 
field  distributions  for  plane  H and  E-wave  diffraction  by  dielectric  cylinder  with  8=4. 
k()r  = 3 are  shown  in  Fig. 2.  where  the  direction  of  incidence  is  indicated  with  arrows. 

The  solution  very  close  to  reference  eigen  function  solution  can  be  obtained  in  wide 
range  of  a here  - in  case  of  E-wave.  The  criterion  of  the  best  choosing  a is  under 
investigation  now. 


Fig.  2 

The  other  features  of  the  procedure  arc  quite  traditional  for  finite  methods:  the 
absorbing  boundary  conditions  are  imposed  upon  the  scattered  field  at  the  mesh  border, 
the  scattering  diagram  is  calculated  with  the  Kirchhoffs  formula  and  from  the 
assumption  of  far  zone  observation  point  location  [3]. 

There  exists  an  opportunity  of  immediate  visualization  of  field  distribution  within  the 
calculation  area,  because  such  a distribution  is  just  the  output  of  wave  equation  solution. 
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For  example,  distribution  of  the 
field  near  and  within  the  illuminated 
in  normal  direction  dielectric  slab 
with  £ = 4 , length 

kQl  = 60  (« 10A,o)  and  thickness 

X0  / 4 = L£  / 2 (half- wavelength 

passing  through  without  reflection) 
is  shown  in  Fig.3.  Wavy  disturbance 
of  the  field  near  the  edges,  regular 
distribution  within  the  slab  resulting  from  interaction  between  going  through  and 
transversal  standing  waves,  can  be  seen. 


Fig.3 


The  backscattering  RCS  pattern  of 
k0l  - 60  («  10A,o)  metal  plate, 

calculated  in  H-wave  case,  is 
presented  in  Fig. 4. 

The  advantages  of  the  above- 
described  procedure  are  its 
universality,  the  simplicity  of 
interface,  the  opportunity  to  introduce 
the  inhomogeneous  media. 


1000o  -g  k0c  Backscattering  RCS  pattern. 

i Metal  plate  of  width  60  («1  Ok)  and  thickness  tt/4. 
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Fig.  4 
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